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Abstract 

Using variational mean-field theory, many-body dissipative effects on the threshold law for quan- 
tum sticking and reflection of neutral and charged particles are examined. For the case of an ohmic 
bosonic bath, we study the effects of the infrared divergence on the probability of sticking and ob- 
tain a non-perturbative expression for the sticking rate. We find that for weak dissipative coupling 
a, the low energy threshold laws for quantum sticking are modified by an infrared singularity in 
the bath. The sticking probability for a neutral particle with incident energy E ^ behaves 
asymptotically as s ~ ; for a charged particle, we obtain s ~ Thus, 

"quantum mirrors" -surfaces that become perfectly refiective to particles with incident energies 
asymptotically approaching zero- can also exist for charged particles. 

PACS numbers: 68.43.Mn, 03.65.Nk, 68.49.Bc, 34.50.Cx 
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Since the very early years of quantum theory, theorists have considered the interaction 
of low-energy atoms and molecules with surfaces In comparison to a classical particle, 
a quantum particle at low energy was predicted to have a reduced probability to adsorb to 
surfaces. The reason is despite the long-range attractive van der Waals interaction between 
a neutral particle and surface, at sufficiently low energies, quantum particles have little 
probability of coming near the surface |2|. 

This effect is named "quantum reflection," and it is a simple result of the wave-like 
nature of low-energy particles moving in a finite-ranged attractive potential. This reduction 
in the particle's probability density near the surface leads to a reduction in the transition 
probability of the particle to a state bound to the surface. In one of the earliest applications of 
quantum perturbation theory, Lennard- Jones and Devonshire concluded that the probability 
of a neutral particle with energy E sticking to the surface should vanish as \/E as — )■ 0. 

In contrast, charged particles do not experience the effects of quantum reflection. Far from 
the surface, charged particles interact with the surface through a Coulomb potential. Due to 
the slow spatial variation of the Coulomb potential, incident particles behave semiclassically. 
As a result, Clougherty and Kohn [2] found that the sticking probability should tend to a 
non- vanishing constant as — )■ 0. 

The seemingly universal scaling law for neutral particles was shown to hold even within a 
non-perturbative model that includes arbitrarily strong quantum fluctuations of the surface 
2|, |3|. This model however was regularized with the use of a low-frequency cutoff. Thus the 
effects of an infrared divergence involving low frequency excitations were not included in the 
analysis. 

In the eighties, experiments went to sub-milliKelvin temperatures to look for this thresh- 
old law scaling in a variety of physical systems without success Q] . Theorists [5| realized that 
the experiments suffered from unwanted interactions from a substrate supporting the target 
of a superfluid helium fllm. By increasing the thickness of the fllm, the next generation of 
experiments produced data consistent with the a/E law, and the controversy subsided. 

In recent years, with dramatic advances in producing and manipulating ultracold atoms, 
there is renewed interest in interactions between low-energy atoms and surfaces. New tech- 
nologies have been proposed that rely on the quantum dynamics of ultracold atoms near 
surfaces; microfabricated devices called "atom chips" would store and manipulate cold atoms 
near surfaces for quantum information processing and precision metrology Our under- 



standing of device performance will depend in part on our understanding of ultracold atom- 
surface interactions. Experiment is now in a position to test detailed theoretical predictions 
on the behavior of low-energy sticking and scattering from surfaces. 

In this Letter, we consider theoretically a non-perturbative model that focusses on the 
effects of low-frequency excitations on quantum reflection and sticking. We follow the mean- 



field variational method introduced 
dynamics of the spin-boson model 



jy Silbey and Harris in their analysis of the quantum 
8|. Using this method we analyze the effects of the 
infrared divergence on the sticking process. Our analysis reveals two distinct scaling regimes 
in the parameter space in analogy with localized and delocalized phases in the spin-boson 
model. In the delocalized regime, an infrared divergence in the bath is cutoff by an energy 
scale that depends on the incident energy of the particle E. As a consequence, we find that 
both the threshold laws for neutral and charged particles are modified by the dissipative 
coupling strength a. As a result of the low frequency fluctuations, the threshold law for 
neutral particles is no longer universal, and the threshold law for charged particles no longer 
precludes perfect reflection at ultralow energies. 

We take for our model a particle coupled to a bath of oscillators 

H = Hp + H, + H, (1) 

where 

Hp = Eclck - Ebb^b, (2) 

iffe = ^Wgajag, (3) 
q 

Hc = -iiclb + b^Ck)gi'^(riujg) (a, - aj) - i^Cfc^fa ^ a (w^) (a, - aj) 

-ib^bgs ^ a (uq) {aq - aj) (4) 
g 

where gi, g2 and g^ are model coupling constants and a (uq) depends on the specific particle- 
excitation coupling mechanism. c|. (c^) creates (annihilates) a particle in the entrance chan- 
nel \k) with energy E; b'^ (b) creates (annihilates) a particle in the bound state \b) with 
energy —Ef,. (a^) creates (annihilates) a boson in the target bath with energy Uq. (We 
use natural units throughout where h = 1.) We work in the regime where E <^ Ef,. We 
neglect the probability of "prompt" inelastic scattering, where bosons are created and the 



particle escapes to infinity with degraded energy, as the phase space available for these pro- 
cesses vanishes as — )■ 0. Thus only the incoming and bound channels are retained for the 
particle. 

We consider a model with ohmic dissipative spectral density. Such a model can be realized 
with a semi-infinite elastic solid where the incident particle couples to the surface strain. 
The spectral density function that characterizes the coupling to the excitation bath is given 
by 

Ji^) =^9l^^ i^q)^i^ -^q) = alp^ (5) 

where p is a frequency-independent constant. 

This model differs in an important way from the model of Ref. ^ where low frequency 
modes were cutoff to prevent an infrared divergence in the rms displacement of the surface 
atom. In this model, low frequency modes are included, and their effects on quantum 
reflection and sticking are the focus of this study. 

We start with the variational approach used by Silbey and Harris jsl for the ohmic 
spin-boson model. A generalized unitary transformation U = e*^ is first performed on the 
Hamiltonian H, with 

S = ib^bx (6) 



and 

x = y'— (cig + aj) (7) 

The variational parameters to be determined are denoted by fq. The unitary transformation 
displaces the oscillators to new equilibrium positions in the presence of the particle bound to 
the surface and leaves the oscillators unshifted when the particle is in the continuum state. 

The transformed Hamiltonian H is given by 



= Hp + Hh + He 
4 



(8) 
(9) 



where 



Hp = Eclck - E^h\ (10) 
He = -iclh ^ giq{aq - a\)e-''' - ib^Ckc''' giq{aq - aj) 



^6 = (12) 



and where gfjg = ^ficx (cug). We define a mean transitional matrix element A 

A^z(e-Y.g,q{aq-al)^ (14) 

where (• • • ) denotes the expectation over the bath modes. 
The Hamiltonian is then separated into the following form 

H^Ho + V (15) 

where V is chosen such that (V) — 0. Hence, we obtain 

Ho = Eclck - E^h'^h - A*clb - Ab'^Ck + ^ c^^aja, (16) 

Q 

V = -clb gig{ag - aj)e-^^ - A* j - b^Ck ^ie^^ giq{ag - aj 

-iclckYg2q{aq - 4) - ib^b^igsq - fq){aq - 4) (17) 

We calculate the ground state energy of Hq in terms of the variational parameters {fq} 
and minimize to obtain the following condition 

which is an implicit equation for /„. For convenience, in the above we have defined 



-at)-A 



e = E + E, = E + E, + y^ '^^''^''" (19) 
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and 



A = v^l^i, (20) 



u = e ""-l, (21) 
fi^ = ^ ^ (22) 



Under the condition A ^ e, Eq. (fT8|) can be simplified to 



/. = (23) 



where 

A2 



z = — 
e 



(24) 



Using Eq. (l23|) . Eqs. (l2T]) . (122|) and (fT9!) can be rewritten in terms of z 

« = (1 + ^)-^?i^e^3V/(i+-) (25) 

^1 = giQzP^c - giQzpz In ^ (26) 
e = E + Eb + glujc - g^pujc : — (27) 

Uc ~r Z 

where coc is the upper cutoff frequency of the bath. According to Eq.(l2^ and the condition 
A e, z Uc must be satisfied. This leads to the following 

^1 ~ gmp^^c (28) 

e^E + E, + glpu, (29) 
Substitution into Eq. (!24|) gives the self-consistent equation for z 

z = K{l + ^)-9lPe'^^p/(^+^J (30) 

where 

X ~ (ampule? .3^. 

E + Eb + gipcoc 

It is straightforward to find the following closed-form expression for z, valid for z Uc 

pK 

z^K{—)^ (32) 

UJc 

where a, the dissipative coupling strength, is given by a = g^p. 
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Depending on the value of a, there are two solutions to the variational parameters fq. 
We see from Eq. [32] that as a —t- 1, 2; —t- 0. Thus, 

,93q « > 1 

/. ^ (33) 



93g 



a < 1 



In the regime where a < 1, we see that the parameter fg for excitations whose frequency 
Uq <^ z vanishes as Uq — ?■ 0. It is this weakening of the coupling to non-adiabatic excitations 
that allows us to extract a finite mean transitional matrix element. In the process, the 
sticking rate is altered from the perturbative result. 

We can now show that the condition A ^ e is satisfied so that our variational solution is 

self-consistent. According to Eq. f l24|) . A/e = \/ zje. For a > 1, z = 0, so A = and A ^ e 

2 

holds true. For a < 1, z ~ fl'i The coupling constant g\ has a dependence on the initial 
energy of the particle E. This can be seen from the transition matrix element 

51, = -^(6,l,|if,|A;,0) (34) 

The amplitude of the initial state in the vicinity of the surface is suppressed by quantum 



reflection. It is a simple consequence of wave mechanics |2| that in the low energy regime, 
g\g ~ ^J~E as — )■ for a neutral particle. For a charged particle, the coupling constant 
behaves as g\g ~ E^l'^ as — )■ 0, as it is not subject to the effects of quantum reflection. 
Thus in either case, the mean-field amplitude A becomes arbitrarily small as E tends to 
zero, while e approaches a non-zero value. Consequently the conditions for our variational 
solution are always satisfied for sufficiently cold particles. 

For A ^ e, the rate of incoming atoms sticking to the surface can be calculated using 



Fermi's golden rule 



i? = 27r^ (6,1, 



k,0 



5(-Eb-E + cOg] (35) 



where denotes a state of one excitation with wave vector q. 
After calculating the relevant matrix elements, the rate becomes 



Jo 



q 




i? = 2vre"^'^^Ui,-^5^:^^ di-E.-E + ujg) (36) 



After some algebra, we find the leading order of the rate R in the incident energy E to be 



R = 2n{-re"glpE, " (37) 

Uc \Eb + abJcJ 

where z given in Eq. |32]is a constant with a power dependence on gi. 

We compare this rate to that obtained by Fermi's golden rule on the untransformed 
Hamiltonian 

i? = 27r ^ I (6, 1, \H,\ k, 0)1' 6 {-E^ - E + Ug) 

= 27rglpE, (38) 

The matrix elements of transformed Hamiltonian are reduced by a Franck-Condon factor 
which gives the non-perturbative rate with an additional dependence on z. 

The coupling constant gi can be expressed in terms of a matrix element of the unperturbed 
states using Eq. [341 We take He to have the general form in coordinate space 

H, = -lU (x) a {ujg) [a, - a\) (39) 

The coupling constant gi is given by 

gi = {k \U\ h) = / <Pl{x)U{x)cp,{x)dx (40) 







(We have assumed the case of normal incidence, however results for the more general case 
follow from decomposing the wave vecto. into .o™aU.d t_e _nt. Q, 
The continuum wave functions have the asymptotic form for a neutral particle 

0fc(x) /c /ii(a;) (41) 

and for a charged particle 

0fc(x) h,{x) (42) 



where k = y/2mE, and hi{x) are functions, independent of E. 

The probability of sticking to the surface s is the sticking rate per surface area per unit 



mcommg particle flux. Hence, s{E) = ^^R. (We use delta-function normalization for 
the continuum wave functions.) From Eq. [37] we conclude that with a < 1 for a neutral 
particle, 

S(E) ~ (43) 



and for a charged particle, 



s{E) ~ C2^"/2(^-"\ ^ ^ (44) 

where Cj are energy-independent constants. 

In summary, we have considered the effects of the infrared singularity resulting from inter- 
action with an ohmic bath on surface sticking. We calculated using a variational mean-field 
method the sticking rate as a function of the incident energy in the low-energy asymptotic 
regime. We have shown that for an ohmic excitation bath the threshold rate for neutral 
particles decreases more rapidly with decreasing energy E, in comparison with the pertur- 
bative rate. We predict a new threshold law for neutral particle surface sticking, where the 
energy dependence depends on the dissipative coupling a. 

The new threshold law is a result of a bosonic orthogonality catastrophe ll| ; the ground 



states of the bath with different particle states are orthogonal. The sticking transition 
amplitude acquires a Franck-Condon factor whose infrared singularity is cutoff by z. As 

n 

with the x-ray absorption edge a new power law results at threshold. The low-frequency 
fluctuations alter the power law to a bath-dependent non- universal exponent. 

For the case of charged particles, we find that dissipative coupling causes the sticking 
probability to vanish as — )■ 0, in contrast to the perturbative result [2]. Thus, "quan- 
tum mirrors" -surfaces that become perfectly reflective to particles with incident energies 
asymptotically approaching zero- can also exist for charged particles. 
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